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Abstract 

de Sitter spacetime and Bunch-Davies vacuum are a solution to the semiclassical Einstein- 
Schrbdinger equations describing the evolution of spacetime geometry and a massive 
scalar quantum field with arbitrary coupling to curvature. The stability of this solution 
is proven by calculating the renormalized energy momentum tensor expectation value 
for small spatially homogeneous deviations from the de Sitter - Bunch-Davies system 
and solving the linearized backreaction problem. A renormalization scheme is developed. 
All momentum integrations are carried out analytically. The general solution is given 
in terms of its Laplace transform. It contains only two artificial instabilities: a constant 
gauge mode and an instability on the Planck time scale lying outside of the scope of our 
semiclassical theory. 
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1 Introduction 



The de Sitter spacetime [1] is of great theoretical as well as cosmological inter- 
est. The former arises due to its high degree of symmetry: with 10 Killing vector 
fields its isometry group 0(4,1) has the same dimension as the Poincare group of 
Minkowski spacetime and therefore the maximum dimension the symmetry group 
of a four dimensional spacetime can have at all. Just this fact makes a lot of cal- 
culations of quantum field theory feasible in the de Sitter spacetime. 
The cosmological interest stems from the exponential growth of the scale factor in 
the spatially flat (k=0) Friedmann-Robertson- Walker parametrization of de Sitter 
spacetime solving some basic problems of the standard cosmology in inflationary 
universe models [2]. 

In absence of a consistent quantum theory of gravity one usually works within a 
semiclassical framework, where the gravitational field is treated as a classical back- 
ground field and only the matter fields are quantized. This is justified as long as 
all relevant inverse time and length scales are small compared to the Planck scale 1 , 
so that quantum gravity effects are expected to be small. 

Since the work of Schwinger it is known that the quantum fluctuations of a charged 
matter field in an electromagnetic background field can lead to the production of 
particle-antiparticle pairs. The same applies to the gravitational background field 
and is known as the Hawking effect. From this observation the conjecture and also 
some claims arose in the literature (see for example [3]- [4]), that in the presence of 
a scalar quantum field like in most inflationary scenarios the de Sitter spacetime 
might be unstable due to particle production and should decay in some sense by 
itself towards a flat spacetime. 

There is no unique observer-independent particle-antiparticle concept in a general 
curved spacetime and different approaches to particle production involving diffe- 
rent approximations led to different answers for this stability question. 

1 in natural units % = c = 1 
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In reference [5] on the contrary this question is adressed in a rather clear and 
reliable manner (see below) based on the energy momentum tensor, an observer- 
independent physical quantity. Unfortunately no sensible results were obtained due 
to technical problems, on which we will comment later. In a further publication [6] 
the main problems were not eliminated. Nevertheless their approach is promising 
and will be adopted as the starting point for the present investigation. 
Within the semiclassical theory the evolution of spacetime geometry is governed 
by the Einstein equations containing the expectation value of the energy momen- 
tum tensor as the source on the right-hand side, whereas the quantum state of the 
scalar field has to obey the Schrodinger equation, which in turn depends on the 
spacetime metric: 

idt\V) = Hg\V) (1) 

The de Sitter spacetime and the Bunch-Davies vacuum, a special state of the quan- 
tum field, are a solution to this semiclassical system of coupled equations. In order 
to investigate the stability of this solution against small fluctuations of the gravi- 
tational field and of the quantum state we will linearize the equations (1) around 
the de Sitter - Bunch-Davies solution. This linearization is the only approximation 
appearing within this work. The linearized Einstein-Schrodinger equations will be 
solved completely and the general solution will be analyzed with respect to insta- 
bilities. 

In the course of this work another publication [7] on the same subject appeared 
also based on reference [5]. We will reach the same conclusions as reference [7] but 
in a more direct way, because in the calculation of the energy momentum tensor we 
will execute the momentum integrations analytically, so that the result is suited for 
a numerical analysis. Furthermore a general coupling of the scalar quantum field 
to curvature will be allowed. 

The paper is organized as follows: section 2 gives a brief review of some important 
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results from reference [8], whereas section 3, the main part of this work, contains 
the calculation of the energy momentum tensor expectation value and the isola- 
tion of its divergencies. The linear stability analysis is performed in section 4. Some 
mathematical tools are collected in the appendix. 

2 Schrodinger picture field theory in k=0 
Friedmann- Robertson- Walker spacetimes 

This section deals with the quantum theory of a free, massive scalar field in a 
spatially flat (k=0) Friedmann-Robertson- Walker (FRW) spacetime, and it is a 
brief summary of some results from reference [8]. 

Since dimensional regularization will be applied later on, we work in d+1 spacetime 
dimensions and on flat <i-dimensional spacelike hyperplanes with coordinates x = 
(xi, . . . ,Xd)- The maximum symmetry of these hyperplanes will greatly simplify 
the following calculations. The k=0 FRW-metric has the form 

ds 2 = dt 2 - a 2 (t) dx-dx = g^dafdx", (2) 

where a(t) is the FRW scale factor. In terms of the Hubble function H(t) := 
a(t)/a(t), a := d t a, one obtains for the Ricci tensor 

R 00 = -d(H + H 2 ) , Rij = a 2 {H + dH 2 )5 i3 , R 0i = (3) 

and for the curvature scalar R = —d{2H + (d + \)H 2 ). 

A scalar field <fi of mass m is supposed to interact only with the classical gravita- 
tional field g^ u and may have an arbitrary coupling £ to the curvature scalar R. Its 
action reads (yjg := a d ) 

S = Jd d+1 x^g l -{g^{d,<t>){d v <t>) - (m 2 -6R)0 2 ). (4) 

The quantum theory is formulated in the Schrodinger picture using a wave func- 
tional to represent the quantum state. This shows very clearly the real time evolu- 
tion character of our analysis. Then the quantum operators are acting in the Fock 
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space on wave functional 

The spacelike derivatives contained in the Hamiltonian can be dealt with by per- 
forming a Fourier transform (dk := d d k / (2tt) c ^ 2 , a*(k) = a(—k)): 

and the Schrodinger equation resp. the Hamiltonian takes the following form: 
id t V[a;t] = HV[a;t] (6) 

For a curved spacetime without an everywhere timelike Killing vector field no 
unique Fock vacuum does exist. Rather there is a whole class of Fock vacua, which 
can all be represented by Gaussian wave functionals. We want our ^[o;;i] to be a 
member of this class. Furthermore we require the quantum state not to break spon- 
taneously the symmetries of the k=0 FRW metric (homogeneity and isotropy of the 
spacelike hyperplanes), which leads to the following wave functional parametrized 
by one function A(k,t) (the inverse Gaussian width): 

tf[ a; t] = N(t)exp(~[d d kA(k,t)a(k)a(-k)-in(t)), (7) 

where N(t) is a real normalization factor, Q(t) a real phase and k :— \k\. 
Substituting (7) in the Schrodinger equation (6) we get the equation of motion for 
A(k,t): 

iA(k,t) = ^jtd-^tj(a- 2 (t)k 2 +m 2 -^R(t)) (8) 

This is Riccati's equation, and by the transformation A(k,t) =: y/g(T(k,t) + 
i^H(t)) it takes the standard form 



if(M) = T 2 (M) + ^H 2 (t) + %H(t) - (a- 2 (t)k 2 +m 2 -ZR(t)), (9) 



which can be converted by T(k,t) =: —id t hiu(k,t) into the linear equation 

u - (f-H 2 + ^H- (a- 2 k 2 +m 2 -ZR))u = . (10) 
4 2 
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From (8) the following equations for A(k,t) can be derived, which are useful for 
the calculation of the energy momentum expectation value: 



-d t 



2 ReA 



1 \mA 



1 \A\ 



a- 2 k 2 +m 2 -£R) 



V9 ReA 
1 



+ 



d t ^fgd t 



1 \ 



2 ReA 



(11) 



(12) 



g 2 ReA y '2 ReA ' 2 Jg 

The energy momentum tensor acting as the source in the Einstein equations is 
defined as variational derivative of the matter action with respect to the metric 
tensor 



SS 



V9 Sg^(x) ' 

Due to the spatial symmetries the expectation value of the corresponding operator 
in the Gaussian state can be written as 



pit) 



\ 



-p{t) 



-pit) 



v 



-pit) 



(13) 



and the explicit calculation leads to the energy density p and pressure p in terms 
of the width A(k,t): 



P 



dk 



( \Mk,t)\ : 



2ReA(k,t)\ g 

dk f\A(k,t)\ 2 
2ReA(k,t) V g 



+ a- 2 k 2 + m 2 + 2i G 00 -2^ — 2 lmA(k, t)) 

V9 / 



+ (--1) a' 2 k 2 -m 2 + 2£a- 2 G n 



+ 4£ (a" 2 k 2 +m 2 -£R) - 4£ - 2£ 2 lmA(k, t)) 



H 



(14) 



G/j,u = Rnu — \ 9iiv R is the Einstein tensor, and dk := dfkj (2n) d . 
Equation (8) shows, that A(k,t) is of the order of k for large k. Hence the energy 
momentum expectation value (14) is quartic divergent and has to be renormalized. 
These ultraviolet divergencies are due to the behaviour of the wave functional 
for field configurations of high energy and momenta (large k) or resp. for small 
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distances and are connected to the local geometry of the underlying spacetime 
manifold. For this reason they should be proportional to local geometric tensors 
which can be absorbed into the gravitational part of the Einstein equations. Thus 
the divergencies can be removed by a renormalization of the physical parameters 
in the Einstein equations (cosmological constant, Newton's constant and addi- 
tional parameters mentioned below). Fortunately the divergencies of the energy 
momentum expectation value can be calculated as a local functional of a gen- 
eral metric tensor by means of the De Witt-Schwinger-Christensen expansion. It 
turns out that in the Einstein equations one has to admit the geometrical tensors 
H^ V) and ^H^, which are the metric variations l/y/g S/Sg^ u of the function- 

al / d m x ^/gR a/3pa R a( 3 P a, I d^x ^fgR 2 and / d^x ^/gR a/3 R a(3 . Their renormalized 
coefficients have to be regarded as additional physical parameters of the theory. 
We will choose them to be zero, since the effects of these terms have already been 
analyzed elsewhere [9]. 

The renormalization scheme consists in a subtraction of the first three divergent 
terms of the De Witt-Schwinger-Christensen series from the expectation value (14): 

(15) 

If |\&) is a state of finite energy density (compared with an adiabatic vacuum as 
will be explained later on), then the divergencies of (\&|T At „| 1 I r ) and (T Mi ,) DSdiv will 
cancel and (T^) ren is finite. It should be noted, that the renormalization scheme 
decides about the physical meaning of the renormalized parameters. 
With dimensional regularization (d — 3 — e) one obtains [10] 

1/11/ m 2 \\ ( -4m 4 4m 2 1 

(T ^ >DSdiv = 16^ U " 2 V 7 + M 4^ J ■ 9vv ~ —l^— 6 ^ 

+ l(^-(^) + (e-i) 2 W^), (16) 

where 7 is the Euler-Mascheroni constant. For our FRW metric (2) the if -tensors 
are explicitly given in appendix B. 
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2.1 de Sitter spacetime and Bunch-Davies vacuum 

As already mentioned in the introduction we use the k=0 FRW parametrization 
of the de Sitter spacetime: a(t) = e Ht)t in (2). H is the Hubble constant and due 
to the maximal symmetry we have = —d Hi g^ v and R = —d(d+ 1) Hi . 
By substituting the conformal time r(t) := e~ Hot /H equation (10) becomes a 
Bessel differential equation: 

{{kr) 2 d 2 kT + krd kT + {krf - v 2 ) u(k,t) = (17) 

where v 2 := d 2 /4 - (m 2 -£i?)/iJ 2 . 

Its general solution is a linear combination of the two Hankel functions H^ and 
H^: 

u(k,t) = B 1 (k) H^\kr) + B 2 (k) Hl 2 \kr) (18) 

A comparison with the adiabatic vacuum in section 3.1 shows that we have to 
require Bi(k) , since our quantum state should have a finite energy density. 
In addition we want the quantum state not to break the de Sitter symmetry. It 
follows that p and p have to be constant over the whole de Sitter manifold. This is 
the case if B 1 and B 2 in (18) are independent of k, as can be seen by substituting 
y := kr in the integrals (14). 

Therefore we end up with B\(k) = and u(k,t) = H^(kr(t)), leading to 

= -H*- l T d H^\kT)H {2 \kT). (19) 

2ReA(k,t) 4 u \ > u \ ! \ ) 

The quantum state specified in this manner is known as the Bunch-Davies vacuum 
[11, 12]. 

Using (11), (12) and (19) it turns out that the integrals (14) involving two Bessel 
functions are of the Weber-Schafheitlin type and can be evaluated analytically. The 
results are given in appendix A. Expansion in e = 3—d yields (ip(x) = T'(x)/T(x)) 

(20) 
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The 1/e-pole can be removed by a renormalization through the subtraction of 
De Witt-Schwinger terms as described in the foregoing section, and we obtain the 
final result: 

+ m 2 (i-\)R+ ±m 2 R-a 2 ) (21) 
with a 2 = R 2 + \{i~\) 2 R 2 ,R=-12 H 2 . 

Due to the de Sitter symmetry this expectation value is proportional to the metric 
tensor and acts just like an effective cosmo logical constant within the Einstein 
equations. Therefore a de Sitter spacetime with its Hubble constant H determined 
by the transzendental equation 

3H 2 + A = %nG NPrcn {H 2 Q ) (22) 

together with the Bunch-Davies vacuum forms a solution of the semiclassical, 
coupled system of equations (1). 

The solutions of (22) for given A have been studied in ref. [13]. Clearly, for every 
given Hq there is a A so that (22) is fulfilled. Thus a de Sitter spacetime of arbitrary 
curvature is possible. 

3 Nearly de Sitter spacetimes 

We are now approaching our main goal: the linear stability analysis of the semi- 
classical solution from the foregoing chapter. Our starting point is that of reference 
[5]: We will consider small fluctuations of the gravitational field and small pertur- 
bations of the matter quantum state. The semiclassical equations are linearized 
around the de Sitter - Bunch-Davies solution. At time t the whole system will be 
given an initial configuration, which slightly deviates from de Sitter spacetime and 
Bunch-Davies vacuum. Then we analyse the time evolution of this deviation. Any 
instabilities would be indicated by (exponentially) growing parts in the general 
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solution for the deviation. 

For simplicity and feasibility only fluctuations of the metric will be considered, 
which do not break its spatial homogeneity and isotropy. This is of course a limita- 
tion, but it has already been shown (for example in reference [14]) that small initial 
anisotropies are damped away by particle production and an automatic isotropiza- 
tion takes place. 

The Gaussian form of the wave functional is not altered by the metric fluctua- 
tions. This is also assumed for its initial deviation. Due to the linearization this 
assumption does not exclude an initial wave functional containing first excitations. 
Moreover the energy density of initial excitations would be subject to the expo- 
nential de Sitter red-shift and have no influence on the long term behaviour of the 
system. 

Firstly we want to compute the change in the Gaussian width of the wave func- 
tional and in the energy momentum components (14) for a given fluctuation of 
the FRW scale factor and initial deviation from the Bunch-Davies vacuum. In the 
sequel quantities related to the unperturbed de Sitter spacetime and Bunch-Davies 
vacuum will get the index 0, whereas a prefix 5 always means the deviation of a 
quantity from its unperturbed value. 

In order to save some ink the Hubble parameter H Q of the unperturbed de Sitter 
spacetime will be set equal to 1 (in addition to H and c). This means that masses 
are measured in units of H . 

Since the results do not depend on the starting time, t = will be used without 
loss of generality. 

Consider now a small deviation of the FRW scale factor from its de Sitter value 



a (t) = e t : 



a(t) 



a (t)(l + I(t)), 



J(t) < 1 



H(t) 



a 



1 + I(t) 



(23) 



a 
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In the following every quantity will be linearized with respect to I(t). The compo- 
nents of the energy momentum tensor are written as 

p(t) = (Too) = Po + Sp(t) 

P(t) = -- d 9 l3 (T t] ) = p + Sp(t), 

where po and po are the unperturbed quantities (21). There are two sources of 
contributions to 5p und 5p: The first one emerges from the explicit appearance of 
the metric in the definition of T^ u , and the other one is due to the dependence of 
the Hamiltonian on the metric leading to a deviation of the wave functional from 
the Bunch-Davies vacuum: 

A(k,t) = A (k,t) + 5A(k,t) 

With the help of (11), (12) and by noting that d t J dk/(2 ReA (k, t)) = we obtain 
from (14): 

Sp = 

+ (m 2 + £ (Roo-R) + d J) 3 t + \ df) J dkS^^^j 

5p = £5(a- 2 R n ) I dk — l — + - r [dkk 2 5(-^—-] 
y s v 'J 2ReA dJ \2ReAJ 

+ {£a- 2 Rn + (i + ai-d))d t + (\-Od 2 ) /^(^) (24) 

Since ReA(k,t) — J y{t) Re T(k, t) we have 

,( 1 \ dl I a-' 1 \ _ 2 (d+2)I _. 2 

<2RS) = -2^A - Ft ' '(site*) = — TB±- a ° F * (25) 

with F t := -^(5^) • 

The Schrodinger equation (9) leads to the following equation of motion for 5T(k, t): 
id t 5T(k,t) -2T (k,t)5T(k,t) = 2 1(t) % 2 k 2 + r(t) (26) 

with 

r(i) := y I(t) + \ i(t) + t6R(t) = \ Of) + i /(t) + &R(t) , 
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From the foregoing section we know that T (k,t) = —i d t \n H( 2 \kr (t)), T (t) = 
a 1 = e - *, so that the general solution of (26) is 

t 

5T(k,t) = -— * [dt'(r(t')+2I(t')a^\t')k 2 )H^ 2 (kr (t')) 

+ (2) 2 1 ' 5T(k,0) (28) 
Using this and (19) we find 

Ffc " % 2(Rer )2 " Ffc 

2 * 

F« := a - d Re(-^//«VoW) j dt\r{t') + 2 1{t') a^\t')k 2 ) H^\kr {t'))^ 

o 

ff> := a^ d Re(^H^ 2 (kT (t))H^ 2 (k)5T(k,0)y (29) 

The initial deviation ST(k,0) from the Bunch-Davies vacuum To(k,0) is part of 
the initial data of the problem. Since a physically meaningful, perturbed initial 
state should have a finite energy density, 5T(k, 0) must have a special high energy 
behaviour: 

+i 

8T(k,0) = 5T {ii \k)-i J2 (ik) n 5T n , (30) 

n=-3 

where 5T^(k) Q faster than A; -3 , and the coefficients 5T n are determined by a 
comparison with the adiabatic vacuum in the next section. Of course ST(k, 0) has 
to be finite for k — > 0. This can be ensured by suitable 5T^ n \k) in the form (30). 

3.1 The adiabatic vacuum 

Although it is not possible to define a unique vacuum state in a general curved 
background spacetime, it is possible to define a state which is vacuos for the high k - 
modes in the limit k — > oo. This can be achieved by using the adiabatic expansion 
of positive frequency for the field modes, which is at the same time an expansion 
in k" 1 and becomes exact in the above limit. We require the width T(k,0) of our 
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quantum state to coincide with the width r ad (/c, 0) of the adiabatic vacuum in the 
limit of large k, so our energy momentum tensor expectation value will have the 
same divergencies as the De Witt-Schwinger one (being also a local expansion). 
This yields a finite renormalized energy momentum density. 

Since the energy momentum expectation value is quartically divergent, the terms 
up to a relative order of k~ A in the Gaussian width r ad (/c, t) ~ k are responsible 
for its divergencies. "In the limit of large fc" means therefore "up to the relative 
order of k~ A in the limit k — » oo" . 

We are now proceeding with the computation of r a d(/c, t) by an adiabatic ex- 
pansion (positive frequency) of the solution of (10). Again the conformal time 
r := — $ dt' /a(t') with dt 2 = a 2 (t)dT 2 is introduced. We substitute u^k^t) =: 
a 1 / 2 (t) Xfc( r ) an d (10) takes the form 

d 2 T Xk(r)+nl(r)xk(r) = (31) 

with Q fc (r) := (k 2 + a 2 m 2 ) 1 ' 2 , m 2 := m 2 - £ C R . 
The adiabatic solutions of positive frequency are 

1 } 

Xk (r) = -===exp(-i d^Wkir')) , (32) 

where the Wk have to obey the following equations: 

1 W" 3 W' 2 
W 2 + - — ^ - - — - Q 2 (W) 
Wk + 2W k 4 W 2 ~ i2fe {66) 

These are solved iteratively order by order (' = d/dr): 

w k = w® + wP + wP+ ■■■ 

= n k 

(2) 3 fi£ _ 1 fi£ 

k s n 3 k 4 n 2 

(4 ) i nj iv) 5 n'i'n' k 13 99 ngng 297 ng 

fc 16 ^ 8 ftf 32 ftj> 32 128 q7 

ft- fV fV ft- A/ 
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Using the above Q k and expanding with respect to A; 1 up to the fourth order 
relative to the leading one we obtain: 

W k = k + j — - - ± (m 4 + 2m 2 (H + 3H 2 ) - £ C (5HR+R)) + O(k^) (34) 

wjf* and the first term in are already of the relative order of k~ 5 and do not 
appear in (34). Using dt/dr = —a(t) the derivatives with respect to r have been 
converted to those with respect to t. Putting things together we get the adiabatic 
width: 

r . jUM) i w . ( H i w k , 

r ad (M) = -i — — = -W k -t{— - - — ) 

Uad(M) a v 2 2 W k > 

! 1 (m 4 + 2m 2 (H + 3H 2 ) - £ C (5HR+R)) + C(A;- 4 ) (35) 



a" 3 A; 3 8 

This has to be compared with our Bunch-Davies width F (k,t). 

Using the asymptotic expansion of the Hankel functions 

ff< 2 )(;) ~ 2; 2^ + 2^3 + fe" + 1 ' 

as well as | — z/q = tjIq/Hq we obtain: 

r (*, «) = „-* - \h + -^f + -±- 2 *f - -L^f K + rhJ) + (0 

(36) 

For the sake of clarity we did not replace if by 1 in this formula. The asymptotic 
expansion (36) coincides with (35) up to the order given in the special case of 
de Sitter spacetime (H = Ro = Ro = 0): T = r ad0 + 0(k~ 4 ). This means that 
Bi(k — > oo) — > was the correct choice in section 2.1. 

For the nearly de Sitter spacetime we have to require 5T(k, 0)+T (k, 0) = r a d(A;, 0) + 
0{k~ A ). After linearizing (35) with respect to the deviation from the de Sitter 
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spacetime we are finally in the position to obtain the coefficients 5T n needed in 
(30): 

ST, = -1(0), 5T = 1/(0), ST.! = -±(rhp(0)-Z c 6R(0)), 
5r_ 2 = ^(Arhll(0) + 2mli(0)-U25R(0)+5R(0))), 
5r„ 3 = -l(3/(0)(77^+6mg)+27T^(/(0)+67(0)-^i2(0)) 

o 

- £ c (68R(p)+58R(p) +5R(0))) . (37) 



3.2 Momentum integrals and isolation of divergencies 

If we insert (25) into equation (24) there are two integrals of the Weber-Schafheitlin 
type involving two Hankel functions, which already appeared in section 2.1. After 
substituting y := kr they are evaluated using (A.l): 



4 2) := fdk i— 

J 2ReA (k, 



t) 



m /2 

- - 7 + 1 + ln47r - ^(§ + i/ ) - 4>(l-vo) + O(e) 



167r 2 Ve 

7(2) ._ (JL afWk 2 , 
J2 - J dk 2ReA (k,t) (38) 

ml 3 , (2 



16tt 2 



- (m 2 -^o) (- - 7 + g + ln4vr - ^(§+i/„) - ^(§-*b) + O(e)) 



The other integrals appearing in (24) are involving four Hankel functions. In terms 
of 



= ^f^/P ] d W l+d - 1 H^(yro(t-t')) H^\y) (39) 

they explicitly read 

2 * 

JdkFjp = Re(-i^ Jdt' (rit')4 4 \t-t') + 2I(t')4 4 \t-t'))^ (40) 
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jd~kao 2 k 2 Fl: 
IdkFi 1 



Re 




o 



Re 



7T 2 1 \ 

i — 5T n i n Ji A) (t) + SVW -terms (41) 




Re 




r 2 (t) £ 5r n i n ji%(t))+5rW .terms. 



n=-3 



The i' -integrations in (40) are convolution integrals. The fact that the j\ ' defined 
above depend only on the difference t—t' means independence of the starting time 
(t = here) and follows from the maximum symmetry of de Sitter spacetime. 
The integrals (39) involving a product of four Hankel resp. Bessel functions were 
not found in the mathematical standard literature. Therefore their evaluation has 
been included as part of this work in appendix A. 

At this point the investigations in references [5] and [6] failed. In order to circumvent 
the integrations of four Hankel functions they carry out a so called "short-time" 
approximation, which consists in a restriction on short times t and the Taylor 
series expansion of j/ 4 ^(£ — i') around t' — t up to the linear order. However, the 
J, (4) (*-f) have a singularity at t' — t (see below) and cannot be expanded around 
this point. Hence this approximation does not lead to correct results even for 
arbitrarily short times. Moreover the divergencies are not coming out correctly, so 
no sensible renormalization is possible. Comparing with our exact results (65) and 
(66) it turns out that even the leading terms of a real short-time approximation 
are missing. 

In reference [7] the momentum integrations are not executed and one is left with 
even more complicated integrals in the final result. 
We show in appendix A Eq. (A. 10) that 





G(t;0,l). 



(42) 
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The dependence of G and on d and v has been suppressed in favour of a shorter 
notation. The definition of G in appendix A explictly reads: 
-r d (t) 



G(t; P ,l) := -TM 
4 sin nu \ 



p fd+l d+l d+l 1 + l 2m 

- 2 4^3^ P,-^ + -; 1+^,1-^, — g — ;T oW 

+ (*) e 4^3 ~2 — ' 2 _iy; 1_Z/ ' 1_2 ^' ~~ 2 ^ T ° 

+ r 2 ^) e- 2 ^ 4 F 3 (^+,- P , ^, ^ + 2„, 1 + ,; 1+u, 1 + 2,, ^+,; r 2 (t))) 

(43) 

Our generalized hypergeometric function 4 F 3 appearing in (43) is defined by an 
infinite series: 



4 F 3 (a 1 ,a 2 ,a 3 ,a 4 ; fa fa fa, z) := 

T(a 1 ) r(a 2 ) r(a 3 ) r(a 4 ) 



4-^3(^1, « 2 , a 3 , a 4 ; /3 2 , 2) 



^2, T(o;i+n) r(a 2 +n) r(a 3 +n) T(o; 4 +n) n 
^0 n! r(A+n)r(^ 2 +n)r(^ 3 +n) * ' ( j 

where p F q is the function usually called generalized hypergeometric function in the 

mathematical literature. 

According to appendix A the integral (39) is convergent for T (t — t') 7^ 1 resp. 
t' 7^ t in the region 4|Rez/|</ + Re<i<3 (which has always a non-zero extension 
if m 2 — £Rq > 0). Equation (42) gives its analytic continuation on the whole complex 
d -plane. 

In order to investigate the convergence behaviour of the series (44), we need an 
asymptotic expansion of its terms for large n. This can be obtained using Stirling's 
series for the gamma-function (see [15, §13.6]): 

na+n) ~ ^ e -^exp(| 'gy + °(^W))~ <«> 

where B m (a) are the Bernoulli polynoms: 

1 1 3 x 

S (x) = 1, = x - - , B 2 (x) = x 2 - x + - , B 3 (ar) = x 3 - -x 2 + - , 

B 4 (x) = x 4 - 2x 3 + x 2 - , £ 5 (a;) = x 5 - ^x 4 + ^x 3 - -x , ... (46) 
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4 4 



Defining 

B m (ai, a 2 , a 3 , a 4 ; Pi,P 2 ,Pz,Pa) ■= ^B m (ai) -^2 B m(Pi) , 

i=i i=i 

° ■= (47) 
1=1 1=1 

we get 

r(ai+n) T(a 2 +n) T(a 3 +n) T(a 4 +n) 



Qn ■ 



r(P!+n) T(p 2 +n) r(p 3 +n) F(l+n) 

m (-)m+iB m+ i(ai,a 2 ,a 3 ,a 4 ;Pi,p 2 ,p 3 ,l) , . 1 



- n ex P^Z. m(m + l)n™ + ^ n M + i/2 



m=l 



M=4 



n .( 1 + ^i + %i + %0 + ftkO +c)( ^ )) (48) 



The polynoms m arise from an expansion of the exponential function and are 
polynoms of Bernoulli polynoms: 

x ggQ ■ ■) ^ . = 0f _ -B3 ^ _ g| £? 2 03 -04 
U "' J ' " 2 ' 2 ' 8 ' 6 ' 3 ' 48 " 12 12 



~ Bo B^B^t B 2 Ba B\ Bk ~ 

B ^-=^h~ + -T^r + ~ 7^ , 0o := 1 ■ (49) 



The series (44) is always convergent for < z < 1 respectively t > in (43). 
Putting in the arguments cci . . . P 3 from (43) we find a = d + l — 3—p and therefore 
< er < 4 for d = 3, / = 0, 2, 4 and p = 0. This means that in the limit z — > 1 
(£ — > in (43)) the leading terms of the series are the ones of large n, which can 
be computed using the expansion (48). 
Introducing the function 

oo ~n 

F ^,s) = Y.- s (50) 

we obtain 



n=l 71 



4 F 3 {ai,...,a 4 ;Pi,...,p 3 ;z) = ^Q n z n 

n 

*= F(z, -a) + Bi F(z, -(7+1) + B 2 F(z, -a+2) 

+ B 3 F(z, -(7+3) + B 4 F(z, -ct+4) + . . . . (51) 
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According to [16, chapter 1.11] the function F obeys the following relations (B m = 
B m (0) Bernoulli numbers): 

oo r> 

F(z,-m) = mli-lnz)-™- 1 -^ -. (l nz y , m = l,2,3,... 

(to + r + 1) r! 

F(z,0) = F(z,l) = ln(l-z) 

1 — z 

F(l,s) = C(s) Riemann's zeta-function . (52) 
Combining the results obtained so far it follows that 

t\t-t') (t-t')- (d+l - 2) ■ (53) 

This means that due to the behaviour of the integrand at the upper limit of inte- 
gration the convolutions in (40) are logarithmically (d — 3, I — 0), quadratically 
(d — 3, I — 2) or quartically (d — 3, I — 4) divergent, if we suppose r(t') and J(t') 
to be smooth functions different from zero. 

The divergent convolutions in (40) represent non-local contributions to the energy 
momentum expectation value (24). However the divergencies emerge directly at 
the upper limit of integration and involve the functions r and / (respectively the 
spacetime metric) only at the point t' = t. Hence the divergencies are again of a 
local nature as claimed in section 2. 

For the purpose of renormalization we have to isolate the divergencies from the 
integrals (40) in terms of 1/s-poles (e = 3—d). Taking into account the behaviour 
(53) of the convolution kernels it turns out that this can be achieved by performing 
l + l integrations by parts. One obtains divergent as well as finite boundary terms 
and finite convolution integrals. At this point the variable p appearing already in 
(43) (being unused up to now) becomes meaningful: According to appendix A our 
G-function satisfies the relations 

Qd tl -(p+l- l -)^G(t-t';p+l,l) = G(t-t';p,l) (54) 
(\d t , - (p + 2- l -))rl{t-t')G{t-t'-p+l,l) = T*(t-t')G(t-t f ;p,l), 
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which can be used to do the integrations by parts in (40). The variable p then 
stands for the number of integrations by parts performed so far. In this way we 
obtain for the integrals appearing in (40): 



±r(t>)G(t-t>;l,o{ * 



Jdt'r(t')G(t-t';0,0) 

l 'l - Jdt'{lf(t')+r(t'))G(t-t';l,0) 
1 o 
t 

Jdt' I (t')G(t-t'; 0,2) = 
o 

\ (V) G(t-t>; 1, 2) - Iff- G(t-t'; 2, 2) + + ^)G(f-f ; 3, 2) 

-M 



j dt'r(t')T*(t-t')G(t-t'; 0,2) = 
o 

■ (r(f) rp(t-t>- 1, 2) - (pQ + r(0) r#?(f-f ; 2, 2) 



2 



+ (^ + ^(t') + 2r(t'))r 2 G(t-t';3,2) 
' ^ + | m + y W + 6 r(f )) r 2 G(t-t'; 3, 2) 



1 rft'J(t') T$(t-t')G(t-t'; 0, 4) 



' 'l(lf)Tp(t-lf; 1,4) - ^^-^2,4) 



+ ^) + 7^ ro ^_^ ;3) 4) _ m + £ /(O +7(0) r 2 G(t-04,4) 



+ {^P- + \ I\t') + m + 3 /(f)) - f; 5, 4)) 



/ // (v) (^) 5 



(55) 
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The convolutions we are left with in (55) are convergent and finite. The same 
holds for the lower boundary terms (proportional to G(t;p, /)) at t > 0. The upper 
boundary terms proportional to G(0;p, I) are divergent of the order a + 1 = d + 
I — 2 — p, because they contain the functions 



4 F 3 (a 1 ,a 2 ,a 3 ,a 4 ; (3i, (3 2 , fa; 1) = E = + . . . 

n=0 n 

= C(-<r) + B x C(-<7+ 1) + B 2 ((-a+2) + B 3 ((-a+3) + 5 4 C(-a+4) 
+ finite terms . (56) 

In the second line above we have already given the analytic continuation of the 
divergent part of the series 4-F 3 (. ..;...; 1) using Riemann's zeta- function. We have 
taken into account as many terms as necessary for the "most divergent" case 
G(0;l,4) wltha = d = 3-e. 

The zeta-function has exactly one simple pole at s = 1 : 

1 oo 

C(s) = — T + 7+£7n(*-l) B , (57) 

S 1 n=l 

where 7 is the Euler-Mascheroni constant (see [16]). This pole becomes an 1/e- 
pole in (56) and represents the divergencies in the typical manner for dimensional 
regularization. In order to separate them from the finite part we rewrite (56) in 
the form 

M 

4 F 3 (ai, ■ ■ ■ , & 5 1) = 4^ 3 («i, • • • , &) + E B m (ai, . . . , fa, 1) C(™ - ^) , (58) 

m=0 

where we have defined the convergent series 4F3 in the following way: 

4 F 3 ( ai ,...,f3 3 ) := (59) 

00 , M s 

Qo(ai, • • • , As) + E (Q„(ai, • • • , As) - E 5 m («i, • • • , As, 1) ^ m 

n=l ^ m=0 ^ 

The number M + 1 of terms to be subtracted from every term of the series is 
determined by a± . . .j3 3 in such a way that the series (59) is just convergent: M 
is the biggest integer less or equal to a + 1, and in our case for d = 3 we have 
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M = l-p + l. 

In (59) d = 3 may already be substituted. The terms of the series behave like 
1/n 2 for large n, hence a truncation of the series at n = N (for example for an 
approximate numerical calculation) will lead to an error of the order of 1/N. 
The only term in (58), for which the regularization d = 3 — e has to be retained 
until renormalization, is the m = M-term in the zeta-function sum containing the 
1/e-pole. 

Now we define a finite function G(p, I) in the same manner as previously G(t;p, I), 
excepted that t — and the 4F3 's in the definition (43) have to be replaced by the 
corresponding 4 F 3 's. 

The functions G(0;p,l) which are needed for the upper boundary terms in (55) 
can then be expressed in terms of G(p, I) and zeta-functions. Using (43), (58), (49), 
(47) und (46) we obtain in the limit e — > 0: 



G(0;1,0) 



G(l,0) + C(l + e) 



G(0;l,2) = G(l,2) + C(-l) + (l + ^cot7rz/)C(0) 



V2V4 / V4sin 2 vrz/ 48 V SV ; 



G(0;2,2) = G(2,2) + C(0)+(i-e(| + ii/oot7ri/))c(l + e) 
G(0;3,2) = G(3,2) + C(l + e) 

G(0;l,4) = G(l,4) + C(-3)+ Q + 3^cot7rz/)c(-2) 




/407 47 2 

£ 1 v + 

V240 96 



v 



• 4 -57z/ 2 /16 
2 sin 2 itv 
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G(0;2,4) = G(2,4) + C(-2) + (3 + 2iz/cot7rz/)C(-l) 

21 3 2 v 2 V , . / 5/1 2 \ 
5 — + 3lI/COt7TI/ )C(0)+ [-TlT-^J 



/49 11 2 3z/ 2 ,21 13 ox 

£ 1/ 5 h« f " ) C0t7TZ/ C(l + £) 

V32 8 4sin 2 7rz/ V 8 6 ' V V ; 



G(0;3,4) = G(3,4) + C(-l)+(^ + ^cot7ri/)c(0) 

/ll 5 2 /83 z/ 2 5 . N \ A/1 

+ V - £ ^ 1" - WCOt 7TZ/ C(l + £) 

V8 2 V48 4sin 2 7rz/ 2 U v ; 

G(0;4,4) = G(4,4) + C(0) + (^3-£(^ + 2z/cot7rz/))c(l + ^) 
G(0;5,4) = G(5,4) + C(l + e) (60) 

The zeta-functions have the following values 

C(l + e) = -+7, a-n) = neN 
e n + 1 

=► C(0) = ~, C(-l) = ~, C(-2) = 0, C(-3) = ^. (61) 

With the help of (25)-(30), (38), (40), (41), (42), (55), (60) and (61) we are in the 
position to specify the rest of the integrals appearing in (24): 
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m = (J + \ + 1 + \ In tt) 

^ol = ^(-^(t) + /(t)(g--^) 
07r^ V 4 v 24 2 sin 7rz/ y 

- 7(0 m*(J + 3 7 - 3 In 2 + ^(§+i/) + \ - i/)) 

+ r(t) 0) - 27(0 2) + 7(0 G fl (2, 2) - + /(*)) G fl (3, 2) 

- r(0) 1, 0) + 27(0) G R (t; 1, 2) - 7(0) 2, 2) 

+ + /(0)) G fl (f ; 3, 2) - / dt' (r(0 + 2r(f)) G*(f-f ; 1, 0) 

o 

+ / ^' + 37(f) + 47(f)) G R (t-t'; 3, 2) 



+ 2 ]T 5T n G R (t;0,n)-^Re j 'dyy 2 H^\y) H^\ya (t)) ST™ (yo (t))) 

n=-3 - o 

(62) 

The abbreviations Gr(p,1) := ReGr(p, /) and Gii(t;p,l) := ReG(t;j9, /) have been 
used. 

Apart from the 1/e-pole in (62) and (63), which will be removed by renormaliza- 
tion, the functions G R (t;p, I) are divergent in the limit t — > 0. Using explicitly the 
asymptotic expansion (51) and the coefficients ST n (37) it turns out that the t— >0- 
divergencies of the lower boundary terms in (55) cancel the t — > -divergencies from 
our resp. 5r(fc, 0) as it has to be. It is for this reason that we need the 5T n - 
terms and the comparison with the adiabatic vacuum in section 3.1. In the same 
way the finiteness of the first and second time derivatives of 5J^ n needed in (24) 
has been checked for t — > 0. 
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<jju> = — -/ (t f — + — z/ 2 + =— i — ) + /(*)( ^ 2 n — J 

2fin g7l2 y kj \\q 16 sin 2 7rz/ / W V 2 8 8sin 2 7rz/^ 

- 6^(0 (A + _ | ^( t ) 

v 48 4 sin 7rz/ y 8 

+ ^ 7(0 m 2 (^ + (m 2 -£i?o)(|j + 7 - In 2 + \ + ± V(§-"))) 

- r(f) 2) + + r(*)) G fl (2, 2) - + \ r(t) + 2r(f)) G R (3, 2) 

+ 27(0 4) - 7(0 G R (2, 4) - + 1 7(0 + 27(f)) G*(4, 4) 

+ + G fl (3, 4) + + \ 7(0 + H 7(0 + 67(0) G«(5, 4) 

+ r(0) rfotf ; 1, 2) - + r(0)) ^(f ; 2, 2) 

+ ( ! T^ + 1 ^ + 2r(0))7#?ii(t; 3, 2) - 27(0) rp R (t; 1,4) + 7(0) rp R (t; 2, 4) 
" + /(O)) rp R {t- 3, 4) + + 3 7(0) + 27(0)) T^ R (t; 4, 4) 

" (^y^ + 1 1(0) + y HO) + 6/(0)) 7#? fl (f ; 5, 4) 

+ J dt' + 3 f (0 + 11 r(0 + 12 r(0) ^(t-t'; 3, 2) 

o 

- / rff + ^ /™(f ) + | 7(0 + 50 7(0 + 48 i(t>)) rp R {t-t>- 5, 4) 
o 

+1 2 00 

-2 £ *r B 7#7 fl (f;0,n+2) - y Re Jdyy A H^\y) H^\ya {t)) 8V^{ya {t))) 

n=-3 o 

(63) 
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3.3 Renormalization 

Our renormalization scheme consists in the subtraction of De Witt-Schwinger terms 
as was explained in section 2. Again we linearize with respect to the deviation from 
the de Sitter - Bunch-Davies system. With the aid of appendix B equation (16) 
leads to: 

= To^6-^ 7+in ^)) 



—4m 



d- 

= i6^U-2^ + ln te))(— K-iW-D' 

+ 4 d(d-3)(l + di + 2(d-2) i) 
9(J ^ ' 

- {£-D 2 d 2 (41 + 4d I + 2{d+l){d-3) (64) 

W = ^(i-i^ + m^^K-Jxi-d)^^/) 

+ — (3-rf)(/ (iv) + 2d 7 + (d 2 +2d-A) I + 2d(d-2) /) 
+ (C-|) 2 ^(4/ (iv) + 8dl+ {6d 2 -4d-6)I + 2d{d-3){d+l)i 

Together with (24), (62) and (63) we finally obtain the components of the renor- 
malized energy momentum tensor expectation value: 
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5p IeQ = 5p-5pT)Sdiv 

= e (sr 00 - sr) 4 2) + (m 2 + td 2 + d dt + \ d 2 t ) (54% + s4t) 

+ 5J Sv + 6J ?L - ^PDSdiv 

1 /in 2 27 
= (— / + m 2 (3 £ c + 1) / + 2&(/ + 3/) + — U + 12 / 

+ l(/ + 3/ + 2/)-36£ 2 / 

~ 2 

+ ^ £ (^oo - SR) (Hl + u) + ^(f-i/) - 1 - mm 2 ) 

+ (m 2 + 9 £ + 3 (£ + 1) $ + | d 2 ) (fc5i?( 7 + 1 + \ In ^) + 8tt 2 5 J^ n ) 

+ ' + 3m 2 / - \ (m 2 + 1) &5i2 - j| - \ (t + 1 + ^ In ^) 

+ 8tt 2 5J 2 ( 1) (65) 

5prcn = 5p- 5pDSdiv 

= ^ (m 2 (A / + /) + ro V + 12 C 2 (/ + 4/) + e c (^- + 3 / + y / + 20 /) 

-^-(/ (iv) + 6/+ll/ + 6/) 
loU ^ ' 

+ ^1^(^11)^(1 + ^) +^(§-1/) " 1 " lnm 2 ) 

+ (£a- 2 ^n + (f - 20 d t + (i-0 9 2 )(^( 7 + 1 + i In ^) + Stt 2 ^!) 
+ 5 / + 3m 2 / - | (m 2 + 1) ^/? - | - J ( 7 + 1 + ^ln ^) 

+ (66) 
According to d = 3 , £ c = £ — |. 

All divergencies have cancelled and the final results (65) and (66) are finite. This 
fact can be regarded a non-trivial check on the calculation. 

Contrary to the other terms the 5jj^ n 's are non-local functionals of r(t) and I(t), 
because they contain convolution integrals as well as the initial data at t — 0. 
The G-functions appearing in the 8 's in (62) and (63) are defined as convergent 
series and are therefore well suited for a numerical computation. 
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4 The linearized Einstein equations 



Restricting on FRW spacetimes the semiclassical Einstein equations 

Rp,v - ^ 9iivR + A-g^ = 8llG N (7) ll /)ren (67) 

are containing two independent components. They read (remember H(t) := a/a = 

H + i(t)) 

3H 2 + A = 8nG N p (68) 
2H + 3H 2 + A = -SnG N p. (69) 

The index "ren" at p and p will be suppressed and the Hubble constant H will be 
explicitly written out in this section. 

If the sources p and p are specified in advance equation (68) (the 00-component of 
(67)) is no dynamical equation of motion but a constraint on the initial data. It 
plays the role of the Poisson equation in electrodynamics. Here instead the sources 
are themselves functionals of the metric and are reacting on its changes. Therefore 
equation (68) is a dynamical equation in the present case. 

Neither the dimensional regularization nor our renormalization scheme are spoiling 
the covariant energy momentum conservation: 

D,(T»°) rcn = p + 3H(p + p) = (70) 

On that account the equations (68) and (69) are not independent: Every solution 
of (68) is also a solution to (69). Therefore, only (68) is considered in the following. 
In its linearized form it reads: 

I(t) 8nG N r . 

6 ^r = -7rt-w;*] ( 71 ) 

Due to the convolutions contained in Sp this is a linear integro-differential equation 
and may be conveniently solved by Laplace transformation. 
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4.1 Laplace transformation 

The Laplace transform of a function f(t) will be denoted by £[/; s] or f(s) 



£[/;*] := f(s) := |^/(*)e- 



•St 



(72) 



Utilizing 



f(s) = sf(s) - /(0) , /<»>(*) = s n f(s) - S -V(0) 



(0), 



dt'f(t')g(t-t');s 



= f(s) ■ g(s) 



(73) 



one is able to compute the Laplace transform Sp(s) of bp (65). For the SJ^ 's 
we need the transforms Gr(s;p, I) and TqGr(s;p, /) of the functions Gn(t;p, I) and 
Tq (t)Gn(t; p, I) , which were defined in (43) in terms of three generalized hypergeo- 
metric series 4-F 3 (. ..;...; Tq (£)) . Their Laplace transforms are again hypergeomet- 
ric series: 



T 2Q 4 F 3 (...;...;r 2 (t)); S 



E 



-2aH t 



Y,Qn(...;...)e- 2H ° tn ;s 



n=0 



Qui-- •;•••) 



(74) 



1 



5-^4 

2if V 2i?o 



S s 

+ a, . . . ; + a + 1, . . . ; 1) 



2#n 



The particular functions 4-F 3 (. Tq (t)) in the 's (62) and (63) are singular 

for t — > 0, so that the corresponding series (74) do not converge. Therefore all 4F3- 

(2) (2) 
series in oJo fin respectively oJ^fin have to be added term by term before summing 

up the series. The complete 5 's are well behaved for t — > (see the comments 

after (62)), and the term by term addition of the series (74) will be convergent. 

The series Gr(s; 1, 0), Gr(s; 3, 2) and Gr(s; 5, 4) appearing in the convolutions are 

convergent by themselves. For large s they are of the order of s _1 (l + In s) . 

The Laplace transform 5p(s) may be cast into the form 



5p(s) 



H o fx f s \ ?r \ ^ W M 



8tt 2 



«r' ' (75) 
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with certain functions 5pj{s/H ) , f n (s/H ) and g(s/H ) . The function g contains 
the 5F(k, 0)- resp. 5T n - and SF^^k) -contributions. 

The explicit calculation shows 2 that the 1(0)- and the 7^(0) -terms drop out in 
(75). Only up to second derivatives of I(t) at time t — are appearing in 5p(s) 
and thereby in (77) (apart from the 5T n , see below). This means that we have to 
specify the initial data 5T(k,0), 1(0), 1(0) and 1(0) in order to fix the solution of 
(71). Regarding 5T(k,0) we only have 5T^ u \k) at our disposal (see (30) and after- 
wards), whereas (37) has to be used for the 5T n . The values of 1(0) and f( lv )(0) 
required in (37) may be obtained from equation (71) and its derivative at time 
t = 0. However, due to the smallness of Newton's constant it is more advantageous 
to specify 1(0) instead of 1(0) and to determine 1(0) from (71) at t — 0. Then it 
can be easily granted that all are small at t — and that the linearization is 
justified. 

We want to save the explicit statement of the lengthy and in the following unessen- 
tial functions f n (s) and g(s), but at least Spj(s) should be given. The expression 
has been left in a somewhat uncompactified shape, because there is presumably 
not much shortening to gain: 



2 This has to prove true since (65) results from (24) and (64). 
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5 P f(s) = ^s 2 + ^(3^)s + 2^(s 3 + 3s 2 ) + ^^ + 12U^c+5)s 



+ ^(s 3 + 3s 2 + 2 S )-36C 2 s 

+ if ^ + 6S ) ^ + + ^(|-) - 1 " ^ | 
+ (^ + 9e + 3(e + |) S +i S 2 )(-6e c ( S 2 + 4 S )( 7 + l + iln^) 



2 

//,7 - ^ ; A. ^ ' "'V ' - ' 2 ~4i?g. 

5 57 z/ 2 m 2 , /5 „ „, ^ 3 ,., . 3 , ,\ 



4 24 2sin 2 7Ti/ iZ^2 
+ (y + f-6^ 2 + 4 5 ))Gfl(l,0) 



2 G R (1, 2) + s G fl (2, 2) - (- + s) G R (3, 2) 



,s 2 



2 j 

- (s + 2) (1. + S - - 6Us 2 + 4s)) Gnis; 1, 0) + + 3s 2 + 4a) 3, 2)) 

+(I* 2 + #+^ 2 ^(f + V^^)) 



/ 1 , m 2 \ 131 49 9 z/ 4 -97z/ 2 /16 

'(^ +1+ 2 ln 4w)-iir-i6'' - 



sin 2 ixv 



/5 17 2 7v 2 x /3 3z/ 2 15 w , 2 , , 

v 2 8 8sm txv' v 8 2sm 7Tf 4 ' 

15m 2 . /3 (ink \ /47 , „ 1 , . 1 ,\\ 

+ (y + \ - GUs 2 + 4s)) 2) + (£ + l) G K (2, 2) 

-(x + i a + 2 ) ^ 2 )) + 2 ^ 4 ) - s d ^ 4 ) + (y + s ) 6 «( 3 ' 4 ) 



4 2 / lv ' V tv ' y uv ' ' V2 

(X + + 2a) G R (A, 4) + (^ + | s 3 n g2 + ^ ^ 4) 



4 2 ' v ' v 8 2 2 

3 2 
+ (^ + 3 s 2 + 11 s + 12) + I - 6Us 2 + 4s)) r^ R ( S ; 3, 2) 



y + I s * + f s3 + 50 s 2 + 48 s ) tPr(*\ 5 > 4 ) (76) 

For a numerical computation of Spf(s) we need the functions G R (p,l) ((43) and 
(59)) and G R (s;p,l) (74). According to the recurrence properties of the Q n (48) 
they can be easily approximated by a direct numerical summation of their cor- 
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responding series up to the N-th term (provided that they are not alternating). 
Since the terms of the series behave like n~ 2 for large n, one encounters an error 
of the order of A r_1 due to the truncation (in case of Gr(s; p,l) , N > s should 
apply). In order that this really comes true, one has to work with a great numerical 
precision. For an example consider the most difficult case G R (1,4): 

g r (ia) ~ E(^ 3 )-^ 3 )) ~ E^" 2 ) 

n n 

Two terms of the order of n 3 calculated independently have to cancel with an 
accuracy of rT 2 in the large-n terms of (59). Hence a relative accuracy of iv~ 5 is 
required for the calculation of these terms, if the summation of the series up to 
n = N is supposed to be sensible. Choosing N = 10000, a relative precision of 
10 -20 is required! 

For that reason we used REAL* 16 variables (29 significant digits). The precision 
approximation for the gamma-function of Lanczos [17] has been extended to a 
relative accuracy < 10~ 25 in the whole complex plane and was used for calculating 
the Q . 

4.2 Stability 

The Laplace transform of equation (71) reads {M Pl = G N 1 ) 
. sT(s) -7(0) H 2 (. i s x f . . 2 , / s n /(»)(0) 1 / s x\ ,„ n 

and is solved by 



fi M_i_ H o fv 2 f ( M J( " )(0) I 1 n( s )\ 
J, , Hp + ^Mf l \^n=QJn{ Ho )-ijnTr + Hq 9Khq)) 



'PI 



Inverting the Laplace transform we get the exact general solution I(t) of the lin- 
earized "backreaction problem": 

+ioo+a 

'« = / 2S r(s)e " (79) 

—ioo+a 
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The contour of integration runs on the right of all poles of the integrand (a > . . .). 
The numerator of I is of the order of s 2 (l + Ins) for large \s\ (on the negative 
real axis between its poles) and the denominator is of the order of s 3 (l + In s) (see 
below). Therefore one can close the contour of integration for t > by a sequence 
of semicircles in the left half complex plane extended by a to the right with their 
centres being located at the origin (compare with the reflected image of figure A.l). 
Their radii have to be choosen in such a way that they do not come too close to 
the poles of the numerator of I on the negative real axis. Then the contributions 
of these semicircles to the integral vanish in the limit of infinite radii. Hence the 
integral (79) is equal to the sum over the residues of all poles of its integrand. 
We are not going to perform the back-transformation explicitly, but investigate the 
question of the existence of any instabilities. Due to the factor e st growing terms 
in I(t) for large t are only possible through poles at s = s with Res > 0. By 
the functions Gr(s;p,1) the numerator of / has infinitely many poles (see (74)), 
but they are all located in the left half complex plane. Thus only the zeros of the 
denominator of / are of interest to us. 
In the following we are seeking solutions of the equation 

6 i - S-M^ (80) 

with Res > . Note that the factor Hq/M pi has to be small compared to 1. This 
is a necessary condition for the applicability of the semiclassical theory. In the 
new inflationary universe scenario for example we have H ~ 10 11 GeV (Mpi = 
10 19 GeV, iJ(today) = 1(T 42 GeV). 

First we will consider two hypothetical situations where instabilities could occur: 
(i) 5p f (s) > In this case a solution of (80) would be 



s „ (3 HI 



H ~ 6tt M 2 Pl 

corresponding to an instability on a large time scale s -1 . However the numerical 
calculations show that f3 is always of the order of TV -1 , which is the error due to 
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the truncation of the hypergeometric series at the iV-th term. This observation 
leads to the conjecture (3 — 0, which can be proven by a simple consideration: 
The constant terms ~ s° in 5pj{s) are coming from contributions to 5p propor- 
tional to I(t) which are present even for a constant I{t). A constant / in the scale 
factor (23) of the metric (2) can be removed by the coordinate transform t' = t, 
x' — (1 + iy^x, for which the 00-component of the energy momentum tensor (21) 
behaves like a scalar and remains unchanged! Hence s = must be a solution of 
(80), but it is a pure gauge mode (coordinate transform). 

With Spj{0) = we have indirectly proven the following two, on this level remark- 
able identities (using (76)): 

2 ^ 2 > = £-d^-f (§ + *- 3h2+ §«t + *> + §«H'>) 

16 16 sin txv 

It is also possible to prove these identities on the level of equation (29), but that 
will be skipped here. 

The numerical result (5 — can be considered a non-trivial check for the numerical 
as well as analytical calculation. 

(ii) 5py(s) — (3s 4 This behaviour would lead to the solution 



v (3 > Wp,^ 



pi 



and thereby to an instability on a short time scale (compared with Hq 1 ) that 
still lies in the semiclassical region. However, again the numerical investigation first 
showed that the s 4 -terms in 5pf cancel each other: 

i((l-12£ c )(G*(l,0) - sGRis; 1,0)) - (G R (3,2) - sG^( S ;3,2)) 

-(1-12^(^(3, 2) - sr§G R (s; 3, 2)) + G R (5, 4) - st§G r (s; 5, 4) 
~ ^(s-^l + lns)) 
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The reason is that fourth derivatives P lv \t) do emerge only in the convolution 
integrals (40) from the singular behaviour of the kernels . Using (53) one sees 
that in fact in (24) the corresponding terms from |cf Jdk 5(2 ReA)' 1 and from 
/ dk 5(a~ 2 k 2 /2 Re A) cancel. 

More easily it follows from the energy momentum conservation 5p + 3H (5p + 5p) = 
, that 5p has to contain one time derivative less than 5p (po = = po + po, see 
(21)). 

Summarizing we did not find an instability but again a non-trivial check on the 
numerical as well as analytical calculation. 

Now it's time to stop making hypothesises and to look at the real behaviour of 
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Figure 1: The function f(s) := 5py(s)/(s 3 + 1) for £ = and certain values of 
m2 := to 2 /Hq. For m2 = 300 a peak appears on the left which is due to numerical 
inaccuracies. Contrary to the rest of the curve it depends strongly on the numerical 
precision iV -1 . 
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8pj{s). The function 8pj{s) / (s 3 + 1) has been plotted in the figures 1 and 2 for two 
values of £ and different values of m 2 for real positive s. The curves are growing 
only logarithmically for large s. 

It has been verified and can also be seen in (76) that the slope of 5pf near s = 
does not reach an excessive large numerical value. Since Hq/M^ is small, a large s 
only can therefore solve equation (80). According to our considerations so far and 
to the numerical investigations 8pf(s) behaves like as 3 (In s + j3) for large s, so that 
(80) has the following solution: 



PI 



(81) 



a' depends on a and j3 without receiving overwhelming large numerical values. 
Hence this solution would lead to an instability on the Planck time scale, which is 
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Figure 2: The function f(s) := 5pf(s)/(s 3 + 1) for £ = 1 and certain values of 
m2 := m 2 /Hi 
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apart from the region of validity of the semiclassical theory. Within our semiclas- 
sical treatment we are not able to conclude for an instability on the Planck scale. 
The investigation of this region remains a subject for a future quantum theory of 
gravity. 

For a constant ratio m 2 / the solution (81) does not depend on H (apart from 
logarithmic terms), so that it would be present even for very small H . Merely 
the existence of our present universe seems therefore to exclude an instability of 
this kind within a complete, applying theory. Probably it is just an artefact of the 
semiclassical treatment, as was argued for example in reference [7]. 
For the practical use of our general semiclassical solution it is always possible to 
avoid this Planck mode by requiring the numerator of (78) to have a zero at the 
same value (81), too. This condition is one constraint in the space of initial data 
which were discussed in section 4.1. The general solution then consists of an infinite 
series of exponentially damped modes due to the poles from (78) with Res < 
and a constant mode from the s = pole. The last one is the only one to sur- 
vive for late times, but it corresponds just to a spatial rescaling of the underlying 
de Sitter spacetime with no influence on physical observables (like its curvature for 
example) . 

5 Summary and conclusions 

We have linearized the semiclassical system of Schrodinger equation and Einstein 
equations (1) around the de Sitter - Bunch-Davies solution in order to investi- 
gate the stability of this solution against small fluctuations of the quantum state 
and small, spatially homogeneous and isotropic perturbations of the k=0 FRW 
de Sitter metric. This linearization in the sense of a stability analysis is the only 
approximation appearing in the present work. 

The condition of finite energy density for the initial quantum state has been dis- 
cussed and the Schrodinger equation was completely solved. The expectation value 
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of the energy momentum tensor has been calculated as a functional of the metric 
perturbation. The momentum integrations have been carried out analytically. A 
procedure for the isolation of divergencies was developed. These were removed by 
a renormalization through the subtraction of De Witt-Schwinger terms. 
The linearized semiclassical Einstein equations became a linear integro-differential 
equation for the metric perturbation. After finding out the initial data which can 
be specified the general solution was obtained in terms of its Laplace transform. 
This Laplace transform has been analyzed analytically as well as numerically using 
the necessary high numerical precision. The general solution contains only two po- 
tential instabilities: a constant mode which corresponds just to a constant spatial 
rescaling of the underlying de Sitter spacetime i.e. a pure gauge, and an instability 
on the Planck time scale which is outside of the scope of a semiclassical theory. 
Thus we have shown that de Sitter spacetime and Bunch-Davies vacuum are stable 
within our semiclassical theory! 

The same result was obtained in reference [7] for the special case of minimal coup- 
ling £ = 0. Since there the momentum integrals are not evaluated explicitly, com- 
plicated estimates were necessary in order to achieve this and there is no numerical 
analysis. 

Our conclusion above is in contradiction with some claims existing in the literature. 
In particular in ref. [4] an instability on the Hubble time scale Hq 1 was found. They 
use indeed a different coordinate system (k=+l FRW), but the coordinate lines 
t = const., on which the spatially homogeneous initial data of the perturbation are 
specified, tend to coincide within the k=0 and k=+l FRW parametrizations for 
late times t — > oo. Therefore the answer to the stability issue should be the same. 
Some criticism on appendix A of ref. [4] is already contained in ref. [7] . Furthermore 
the diverging perturbation a{rj) of the conformal factor found in ref. [4, section 5] 
is very near to a pure gauge mode. It can be transformed into our I(t) remaining 
small for all times and contains no physical divergence. 

Finally a few possibilities for a continuation of this work should be mentioned: 
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Fermi fields could be included, but presumably they wouldn't alter the stability 
argumentation. The investigation of the minimally coupled massless case would be 
interesting, because the renormalized two point function in de Sitter spacetime (see 
refs. [12, 18, 19]) as well as our result (76) are containing a logarithmic infrared 
divergence (the linear infrared divergencies of the ip- and G -functions cancel). At 
last one could try to access the region of quantum gravity via the Wheeler-De Witt 
equation. 

I would like to thank W. Buchmiiller for many helpful discussions in the course of 
this work and D. Litim for many valuable comments on the manuscript. 
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A Integrals of products of Hankel functions 



The definite integral of a product of two Hankel functions can be evaluated using 
the Weber-Schafheitlin integral [20]. One obtains: 

OO 

Jd x ^ H m( x )HP( x) = ,-J r (-£)r(^)«™r(±±* + i,)r(i±*-i,) 



(A.l) 

Jdx^\x) = -^-f e -^-f)i^r(l±^ + ,)r(l±^-,) 

(A.2) 

Jdxx^Hfyx) = ^-§e^-t)Iffl|r(^ + ,)r(i±^-,) 

(A.3) 

Some functional relations for the gamma-function (duplication and supplement) 
have been applied. The integrals (A.1)-(A.3) are absolutely convergent within the 
region 2|Rez/| — 1 < Re A < 0. The left inequality stems from the behaviour of 
the integrand for small x and secures infrared convergence. From the asymptotic 
behaviour of the Hankel functions the right inequality follows, which stands for 
ultraviolet convergence. The integrals (A.2) and (A.3) are convergent even for 
Re A < 1 due to the oscillating behaviour of the Hankel functions. 
In section 3 we need the following integral of a product of four Hankel functions: 

oo 

1 := J dxx^H^ax) H^\x) , < a < 1 (A.4) 
o 

Such an integral of four Bessel functions was not found in the mathematical stan- 
dard literature. Hence it will be evaluated explicitly now. The integrand behaves 
like x A_1 ~ 4 l Rci/ l for x — > and like x x ~ 3 for x — > oo, so that the integral (A.4) is 
absolutely convergent in the region 4|Rez/| < Re A < 2. If a ^ 1 it is convergent 
even for Re A < 3 due to the oscillating behaviour of the integrant. We suppose 
that we are within this region since then all integrals and series appearing in the 
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following will be convergent, too. 

For the square of Hjp (x) we use an integral representation: 

-\-ioc — 0£ 

712 J 2m 1 ( ^ — s) 

—ioo-a L 

V x E R + and \Rev\ < a < | . Integrals of this type are called Mellin-Barnes 
integrals. In order to prove (A. 5) we close the contour of integration by a sequence 
of semicircles in the right half complex plane extended by a to the left with their 
centres in the origin according to figure A.l. The radii of these semicircles have 
to be choosen in such a way that they do not come too close to the poles of the 
three gamma-functions in the numerator. Using functional relations and Stirling's 
asymptotic expansion for the gamma-function one shows similar to [15, §14.5] that 
the contributions of the semicircles to the integral are vanishing in the limit of an 



Im s 




Figure A.l: The contour of integration in the complex s -plane is being closed by 
a sequence of semicircles, whose radii tend to infinity. The poles of the integrand 
in (A. 5) are marked on the real axis. 
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infinite radius. Thus the theorem of residues can be used to evaluate the integral 
(A. 5). The gamma-function T(s) has simple poles at all negative integers s = —n, 
n E N with the residues (— ) n /n\ . The integral (A. 5) becomes the sum of the 
residues of all poles of the three gamma-functions in the numerator of the integrand: 



7T2 n=0 n\ 



In — inn 



_|_ x 2(n-v) e -iir(n-v) 



_|_ x "2(n+v) e -in(n+v) 





—n + v) 


r(- 


-n — 








— n 


) 




r( 


— n + v) 
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n- 


- n 4- 
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—n — v) 


r(- 


■n — 


2v) 




n\- 


- n — 







Via functional relations for the gamma-function and the power series representation 
[20, §5.41] of the product of two Bessel functions J v the proof of (A. 5) can be 
completed: 

- 1 ( n_-«n^ (-rr(l + 2n)(x/2) 2 « 



* sin 2 z/vrV ±^ 



+ E 



^ n! r(l - i/ + n) r(l + i/ + ra) T(l + n) 
(-)" T(l - 2i/ + 2n) (x/2) 2n - 2u 



^ n\ T(l-u + n) T(l - 2v + n) T(l - v + n) 
2inu ^ (-rT(l + 2u + 2n)(x/2)^ \ 
^ n\ T(l + i/ + n) T(l + 2i/ + n) T{1 + i/ + n)J 

-1 



sin 2 z/7T 



-2 e"**" ■/„./_„ + J_ 2 „ + e~ 2 ^ J, 2 ) q.e.d. 



With the integral representation (A. 5) our integral (A. 4) reads after interchanging 
the order of integration: 



+ioo— a ,, , . . \°° 



2 /■ <fc o^_ i7rs r(- s ) r(- a - 1/) r(- a + *) 



3 



e / a e 



7i2 J 2m r(| — s) 

-ioo-a V2 ' 

Choosing a such that Re A — 2 < 2a < Re A — 2 \Reu\ we put in the Weber- 
Schafheitlin integral (A. 3) and obtain: 

+ico-a 

7T 3 V ; J 2ni K ' 

— ioo— a 

T(-s) T(-s - v) T(-s + v) T(| + s) T(f + s + v) £(|_+£ - tj 

r(i- s ) r(Y + «) 

(A.6) 



Jdxx x ' 1+2s H^\x) 
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For a < 1 we can close the contour of integration again by a sequence of extended 
semicircles in the right half complex plane with their radii tending to infinity (figure 
A.l). From the theorem of residues we find: 
-7- _ H) A f^r(l + n)( „ 2n T(l-u + n)r(l + u + n)T(l 

n=0 



n) 



7T 2 sin 2 z/7T ^ n! V r(l + z/ + n) T(l - v + n) r(^±± + n) 



a 2„-2, e 2^ r(f-z/ + n)r(f-2z/ + n)r(l-z/ + n) 
r(l - 1/ + n) T(l - 2v + n) T(^±i - v + n) 
, a 2n +2 , e -2^ T{± + v + n)T{± + 2v + n)T{\ + v + n) 
r(l + 1/ + n) T(l + 2z/ + n) r(^±i + 1/ + n) 

/ ~ ^ /A A A 1 1 A + f 2 x 

{-2 4 F 3 (-,--v,-+v,-;l+v,l-v,—;a J 



7T 2 Sin 2 Z/7T 



+ a" 2 V™ 4 F 3 - - 2v, ±, \ - v- 1 - 1/, 1 - 21/, ^ - 1/; a 2 

+ a^e" 2 ** 4 F 3 (^+^ \\ + ^\ + ^ 1+^ 1+2*, ^ + ^; « 2 )) (A.7) 

We have applied functional relations and introduced our generalized hypergeomet- 
ric function 

^ , ^ r(a± + n) ■ • ■ r(a v + n) „ , . „, 

v F q (a u ...,a p ;l3 l ft;,) := E^p^' <A*) 

r(ai) • • • r(a!p) 

where p F g is the usual generalized hypergeometric function. 

Note that (A. 6) is an integral of the Mellin-Barnes type, too, and is therefore a 
special case of the general notion of Meier's G-function [16]: 

1 




, v , —v 

In section 3 we have A = d + I and a = e~* , and it turns out to be reasonable to 

introduce the following function G: 

A -e~ td ( - ( d+l d+l d+l 1 d+l+1 _ 2t \ 

G(t;p,l) := s -2 4 F 3 p, 1/, Yv,-\\+v,\-v, ;e 

V ; 4sin 2 7rz/V V 2 F ' 2 '2 '2' ' 2 ' / 

2i/t 2iriv 7^ fd+Z d+Z „ d+l 1 d+Z + 1 _ 2 A 

+ e 2ut e 2m \F 3 ( — v-p, — 2u, — , - - v- 1 - v, 1 - 2v, — u; e 21 ) 

+ e" 2 ^ e - 2m \F 3 (— + u-p, — , — + 2^, -+1/; l+i/, 1 + 2^, ^ — + e~ 2t ) J 

(A.9) 
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Then our integral (A. 4) reads just 



e- td l = -^-(-i) d+l G(t;0,l). (A.10) 

With the auxiliary variable p it is possible to express the time derivative of a 
G-function in terms of G-functions: 

(ld t +p + l- l -^G(t;p+l,l) = -G(t;p,l) {All) 

This relation can simply be proven by substituting the definitions (A. 9) and (A. 8). 

B Geometrical tensors for the k=0 
FRW-metric 

The k=0 FRW-metric 

5-00 = 1, goi = 0, gij = -a 2 (t)5ij (B.l) 
leads to the Christoffel symbols 

r°. = Ha^ij, r} = r 0j = H5), otherwise, H = a/a. (B.2) 
From them we obtain the Riemann tensor 

p£t f) p/i f) pM _i_ p^ pA pM pA 

-n- vpa ■— u p L ucr U (? L up ' L p\ L uv 1 ctA 1 pv • 

R°iOj = —R°ijo = a 2 (H + H 2 )5ij, R l ooj = —R\jo = (H + H 2 )Sj, 

R\ji = a 2 H 2 (5)5 k i-5\5 kj ) , otherwise, (B.3) 

as well as the Ricci tensor and curvature scalar 

Rpv '■= R X p\v '■ Roo = —d(H + H 2 ), Rij = (H + d H 2 ) a 2 5ij , R 0i = 
R := g^Rpu = -d{2H +{d+l)H 2 ). (B.4) 
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The if -tensors read: 

5 



:= ^ — Jd m x^\R a ^R a p pa 



— ~-^gpuR al3pa Ra/3pa + 2Rpaf3pRu afSP + AUR^ — %R;iiu ~ AR^ a R° 'v 

+ AR a/3 R aiJ ,p v 
(1) H,v := -^=^ fd m xJ\g~\R 2 



i Sgi- 

~-^9pvR 2 + ^RRpv — 2R;pv + ^9pv U R 



(2) ^ == 



101 

- -\gp,R at5 R a p + c 2R^R av -2R^. ua + nR^ + ^g^uR (B.5) 

Their non-zero components are: 

H 00 = d(-4HH + 2H 2 -4dH 2 H + (3-d)H 4 ) 

Hij = a%(4H + 8dHH + 6dH 2 + 4(d 2 + d-3)H 2 H + d(d-3)H 4 ) 

WHoo = d 2 (-AHH + 2H 2 -AdH 2 H - \(d+l)(d-3)H A ) 

w Hij = a 2 5ij d(4H + 8dHH + QdH 2 + (Qd 2 —4d—Q)H 2 H + \d(d+l)(d—3)H 4 ) 

( 2 )# 00 = d{-{d+l)HH +\(d+l)H 2 - d(d+l)H 2 H + \d(3-d)H A ) 

WHij = a%((d+l)H + 2d(d+l)HH+ \d(d+l)H 2 + d(d 2 +3d-Q)H 2 H 

+ ±d 2 (d-3)H 4 ) (B.6) 

The Weyl tensor vanishes for the conformal flat FRW-metric, so that the relation 
d(d—l) H^ u + 2 — Ad^H^ = results. In d+1 = 4 spacetime dimensions 

the Euler number n = jd*x \g\ 1/2 (R^ P aR pupa + R 2 - AR^R^) is a topological 
invariant, hence its variational derivative H^ v + ^H^ v — 4 ^H^ v vanishes. The 
result (B.6) has been checked explicitly against these two identities. 
For the de Sitter spacetime we have H — and all if-tensors are vanishing for 
d = 3. 
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